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Abstract 

We demonstrate one of the possible ways to compensate the large negative quantum-topological 
QCD contribution to the vacuum energy density of the Universe by means of a positive constant 
contribution from a cosmological Yang-Mills field. An important role of the exact particular solu- 
tion for the Yang-Mills field corresponding to the finite-time instantons is discussed. 
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I. INTRODUCTION 



Current accelerated expansion of the Universe is commonly attributed to the existence of 
the so-called Dark Energy which is confirmed in many cosmological observations so far, e.g. 
in studies of the type la Supernovae [l|, cosmic microwave background anisotropies [1], large 
scale structure jsl etc. The Standard Cosmological Model is based on the time-independent 
Dark Energy approximation called the cosmological constant, or A-term, approximation 
which agrees well with current observational data. However, the problem of theoretical 
interpretation and prediction of fundamental properties of the Dark Energy (or the A-term) 
remains one of the major unsolved problem of Theoretical Physics For a comprehensive 
overview of existing theoretical models and interpretations of the Dark Energy, see e.g. 
Refs. jsj-lsj and references therein. 

One of the traditional interpretations of the A-term is by means of the vacuum energy 
density satisfying the equation of state Pa = — Pa with vacuum pressure Pa and energy 
density pa- However, individual vacuum condensates known from particle physics e.g. those 
which are responsible for the chiral and gauge symmetries breaking in the Standard Model, 
contribute to the vacuum energy of the Universe individually exceeding the observable value 
of the A-term density Ae^p = (3.0 ± 0.7) x lO-^^ MeV^ by many orders of magnitude in 
absolute value j^. This situation, which sometimes referred to as the "Vacuum Catastrophe" 
in the literature, requires extra hypotheses about (partial or complete) compensation of 
vacuum condensates of different types to the net vacuum energy density of the Universe (see 



e.g. Ref. llQl). A dynamical mechanism for such huge cancelations and corresponding major 
fine-tuning of vacuum parameters is yet not known and is a subject of ongoing intensive 
studies (for a review on this topic, see e.g. Ref. jsf and references therein). 

Within the general problem of vacuum condensates cancelation, the QCD vacuum contri- 
bution has a special status. Various existing cancelation mechanisms refer essentially to an 
unknown high-scale physics beyond the Standard Model e.g. to Supersymmetry How- 
ever, they cannot be applied for a compensation of the specifically non-perturbative and 
low-energy QCD contribution. In this paper, we focus primarily on elimination of this most 
"difficult" part of the vacuum energy of the Universe. 

In the framework of the popular instanton liquid models the topological (or instan- 
ton) modes of the QCD vacuum (which sometimes referred to as the quark-gluon condensate) 
are given essentially by the strong non-perturbative fluctuations of the gluon and light sea 
quark fields which are induced in processes of quantum tunneling of the gluon vacuum be- 
tween topologically different classical states. The topological instanton-type contribution 



(top) to the energy density of the QCD vacuum Eyac is one of its main characteristics [12 



and can be written as follows (see also Ref. 13| ) 

+ (0| : rrisss : \0)] ~ -(5 ± 1) x 10^ MeV^ , (1.1) 



Svacitop) = - — (0| : ^P;,(x)Pr(a;) : |0) + 7 (0| : m^uu : |0) + (0| : mjd : |0) 



which is composed of gluon and light sea u, d, s quark contributions. Clearly, other contri- 
butions of a different physical nature should compensate the topological QCD contribution 
(11. ip to the vacuum energy of the Universe since its value by far is not compatible with the 
cosmological observations and data on the A-term value [2|. This issue triggers the search 
for possible cancelation mechanisms, and one such mechanism will be discussed further in 
this paper. 
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II. CLASSICAL EVOLUTION OF THE COSMOLOGICAL YANG-MILLS FIELDS 



Consider one of the possible ways to eliminate the microscopic QCD vacuum contribution 
(11. ip to the vacuum energy density of the Universe introducing the hypothesis about the 
existence of the cosmological macroscopic Yang-Mills fields in early Universe. 

Cosmological solutions for classical Yang-Mills fields have a long history referring back to 
the late seventies, when there was an active search for solutions to the Einstein- Yang-Mills 
field equations Later, the role of non-Abelian gauge fields in the early Universe evolution 
has been intensively studied in many different aspects, in particular, in the context of the 
Dark Energy [l5| and non-Abelian fields driven inflation without a presence of a scalar field 
( "gauge-flation" ) [l6|. Practically, there are no any physical arguments which could forbid 
the existence of the homogeneous non-Abelian gauge field with unbroken SU{N) symmetry 
at cosmological scales with an isotropic energy- momentum tensor possibly originating 
from the inflationary stage of the Universe evolution [16] . 

Let us now assume that a chromodynamical (gluon) field with unbroken color SU{3)c 
symmetry exists as a real physical object filling up the early Universe, and the subject 
of our further discussion concerns possible physical states of this field and their real-time 
dynamics. For simplicity, we work in the fiat Friedmann Universe with conformal metric 
dfiu = o?{.il)9nv{M)^ where g^u(M) is the Minkowski metric. The Einstein equations with 



energy-momentum tensor of classical Yang-Mills fields are [14 




:0, F;^ = d,Al-d,Al + r"^AlAl. (2.1) 



-9 J \/-9 

This system is written in the most trivial form without taking into account interactions of 
the macroscopic Yang-Mills field with the physical vacuum (no vacuum polarisation effects 
are included here) and other forms of matter (i.e. e = 0). Here and below, raising and 
lowering Lorenz indices are done by the Minkowski metric 9^u(m) as usual. 

Since initial conditions in the early Universe are quite arbitrary, it is meaningful to 



start with the study of spatially- homogeneous and isotropic modes of the gluon field [ITJ . A 
specific feature of such modes concerns their distinct topological structure where the isotopic 
and spatial indices are mixed up. In the case of Hamiltonian gauge = and homogeneous 
and isotropic 3-space we have the following simple structure of these modes: 

^a^[ ^tm. ^, a = 1,2, 3 
^ \o, i = 1,2,3; a > 3, ^ ' ' 

with a single non-trivial time-dependent degree of freedom A{ri) to be studied in what 
follows. In this case, the classical Yang-Mills equations (12. ip read 

= [A" + A') , A" + 2^3 = , (2.3) 

and thus completely determine the conformal time evolution of the spatially-homogeneous 
and isotropic Yang-Mills field. The second equation in Eq. (12. 3p can be exactly integrated, 
and its general solution implicitly corresponds to non-linear oscillations, i.e. 
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with C, Aq being integration constants. Numerical solution of Eq. (12. 3p for the gluon field 
potential with initial condition A'{0) = and an arbitrary amplitude Aq = C to a good 
accuracy can be approximated by 

A{r]) Ao cos (^^Ao^y (2.5) 

An essentially non-linear character of oscillations of the classical YM field is thus emerged 
in explicit dependence of their amplitude on frequency. According to Eqs. (12.31) and (12.41) . 
the spatially-homogeneous classical YM field in the isotropic Universe behaves as an ultra- 
relativistic medium with energy density ^ym ~ and equation of state pym = ^ym/3 
a. 



III. ROLE OF THE VACUUM POLARISATION 



Can a classical Yang-Mills field be a component of the cosmological medium in the 
radiation-dominated Universe? A simple analysis have shown that the classical spatially- 
homogeneous Yang-Mills field cannot exist in the early Universe since the classical Yang-Mills 
equations (12.11) are not form-invariant and unstable with respect to radiative corrections. 
Such an instability emerges due to the fact that there is no any threshold for vacuum po- 
larization of a massless non-linear gauge field, i.e. any infinitesimally small external field 
is capable of reconstruction of the classical Yang-Mills vacuum [l8|. Due to non-linearity 
of initial operator Yang-Mills equations the vacuum polarization of the massless quantum 
gluon field by its classical component leads to a modification of classical equations. In prac- 
tice, we deal with the Savvidy equations for the Savvidy vacuum fluctuations 18| and look 
for their spatially-homogeneous modes. 

Next, let us analyze the Yang-Mills equations incorporating the vacuum polarisation 
effects. The Lagrangian of the gluon field taking into account the vacuum polarisation in 



the one-loop approximation has the following form [18 



4^YM 



2 



11 



1287r2 



■^^In 
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J 



^^QCD 



J 



^^QCD 



J 



(3.1) 



Here, the numerical parameter C, is not fixed and refiects an ambiguity in normalisation of 
the corresponding gauge/Lorentz invariant J. Such a Lagrangian leads to a modified system 
of equations for gravitational and Yang-Mills fields in the isotropic Universe with vacuum 
polarisation effects incorporated, namely. 



rpu,maX, _|_ \xy I ^ 



o-A 



^ah 



duV^9 - r'^'Ai 




pa puX 



(3.2) 
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where e ~ 2.71 is the base of the natural logarithm; Aqcd is the QCD energy scale; J"^'™*^* = 
{e + p)u^u'^ — d^j) is the energy-momentum tensor of all components of the cosmological 
medium except for the macroscopic Yang-Mills field; A = Ainst + Acosm + • • • is the total 
contribution to the vacuum energy density which consists of the non-perturbative spatially- 
inhomogeneous (topological) quantum fluctuations of the gluon and quark fields (quark- 
gluon condensate) of an instanton nature (11 .ip . Ajnst = ^vac{top) — "(5 ± 1) x lO^MeV^, 
an uncompensated contribution from the observable cosmological A-term, Acosm, and dots 
represent all other perturbative vacua contributions. The A-term value, Acosm, could have a 
different nature, other than topological non-perturbative one in QCD or perturbative ones 
in high-energy particle physics, so we explicitly separated it from the rest. Form now on, 
we implicitly assume that perturbative components of the net vacuum energy density from 
all other microscopic vacuum condensates in particle physics are compensated elsewhere at 
high energy scales and do not enter the vacuum energy density of the Universe, so A = 

Ajnst ~l~ Acosm- 

The system of equations (13. 2p is written in the most general form including all forms of 
matter, as well as the uncompensated quark-gluon condensate contribution Ajnst and the 
observable cosmological A-term Acosm- The components of the energy-momentum tensor for 
the homogeneous and isotropic modes specified in Eq. (12. 2p have the following generic form: 



,0,tot _ „0'™at T 33 1 



T^'- = To"' +A + 



647r2 



f,p\ A''' — 4^ 
{A" + A') In ^"'^^ ""J + A'' - A 



IS, tot ^ ^/3,mat , T r/3 , _}]_]_X0 



a^(eAQCD)^ 



(3.3) 



In fiat and isotropic Universe, trace of the Einstein equations and the equation of motion of 
the macroscopic gluon field read, respectively, 

= . - 3p + 4A + , T-^^ = -1^1 {A- - A') , (3.4) 



^ ( A' In ^ ) + 2A^ In ^ = 0. (3.5) 

^^ V «'(eAQCD)^ ; a4(eAQCD)^ ^ ^ 



It is straightforward to show that the (0, 0) Einstein equation 



3 a'2 - 33 1 
7 = £ + A + 



(A'2 + A^) In ^ + - A' 

^ ^ a4(eAQCD)^ 



(3.6) 



is the exact first integral of the system of equations (13. 4p and (13. 5p . while the exact first 
integral of second equation (13.50 is 

6e{A'^ - A^) _ 

a^(eAQCD)^ ^^-^^ 
The latter leads to a considerable simplification of the energy-momentum tensor, namely, 

„0.tot _ ^o.mat T 33 (^Aqcd)^ 

^° +^+ 647r2 6e ' 

rr' = rr^' + (a + ^}^^) € - (3-8) 
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Now we can observe an interesting possibility to eliminate the microscopic negative QCD 
contribution to the vacuum energy, Ainst, by means of the constant positive contribution 
from the spatially-homogeneous mode of macroscopic gluon field. A small non-compensated 
remnant - the observable A-term, Acosm ^ can, in principle, have a different nature which will 
be discussed in our forthcoming publication. The corresponding condition for the Ajnst — 
—265^ Me compensation 

^^^^^f^ + Ainst = 0, AQCD^ 280MeV, (3.9) 
647?^ be 

however, is not fully automatic; it is satisfied for a certain value of the normalisation pa- 
rameter only, ^ ~ 4, which should be constrained in a complete theory of the QCD vac- 
uum. Therefore, in principle, one succeeds to eliminate the huge negative contribution from 
spatially-inhomogeneous non-perturbative quantum fluctuations of the gluon field by means 
of a positive contribution from fluctuations of spatially-homogeneous macroscopic gluon field. 
This is achieved by fixing the remaining freedom in normalization of the Yang-Mills invariant 
J in the Lagrangian f l3.ip . As we will see below, both mutually compensating contributions 
to the vacuum energy density of the Universe have a common instanton nature. 



IV. COSMOLOGICAL EVOLUTION OF FINITE-TIME INSTANTONS 

Together with the compensation condition (13. 9p and the first integrals (13. 6p and (13. 7p . 
the resulting system of equations (13.40 and (13. 5p is dramatically reduced to the following 
simple form: 

3 a'^ 

-=e + Acosm, (4.1) 

X a* 

A" -A' = a^Mo^^ . (4.2) 
6e 

Notice that under the exact cancelation condition (13.91) the cosmological (macroscopic) evo- 
lution of the Friedmann Universe given by the scale factor a = a{ri) is now completely 
decoupled from the microscopic evolution of the gluon field A = A{t]). The physical time 
scale for the cosmological evolution is of the order of the Universe age tcosm ~ ^/H (in terms 
of the Hubble parameter H), while the typical time scale for the Yang-Mills field evolution 
is of the order of the hadronisation time thadr ~ 1/Aqcd- So at present epoch the right hand 
side of Eq. (14. 2 p can be taken to be constant in time to a good accuracy, or more precisely, 
given by a classical solution of the Friedmann equation (14. ip . In practice, this means that 
the dynamical cancelation under the condition (13.90 and, hence, the decoupling of the QCD 
vacuum fluctuations from the hot cosmological plasma have effectively happened at the end 
of the hadronisation epoch in the early Universe evolution. 

For convenience, let us rewrite the Yang-Mills equation (14. 2 p in terms of dimensionless 
time and gauge field as follows 

§r-A^^l, A^A^fl^-^, (4.3) 

dv J ^Aqcd 2Aqcd (6e)V4 

for ^ — 4. For simplicity, we have chosen the initial values of the Yang-Mills field and the 
scale factor as follows: 

A{fio = 0)=Ao = 0, a{f,o = 0) = ao = 1 . (4.4) 
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We can see now that, indeed, the time scale of the Yang-Mills field fiuctuations is essentially 
microscopic and corresponds to the Aqcd energy scale. The equation (14.31) can then be 
easily integrated, and its general solution can be written in the following implicit form: 

^ dA _ 

h (4.5) 



where Aq is an integration constant. Notably, the analytical solution (14. 5 p taking into 
account the QCD vacuum polarisation, in fact, differs from the classical Yang-Mills solution 
(12. 4p by sign in front of under the square root only, having though a significant effect 
on its time dependence. Moreover, since the solution (14. 5 p was obtained under the exact 
cancelation condition (13. 9p . it corresponds to the minimal energy of the QCD system in the 
ground state of the Universe and, hence, is physically preferable. 

For the initial conditions given by Eq. (14.40 (independently on a particular do value), the 
solution for Aifj) (14.50 obeys the following properties: 

• Symmetry: A{—fi) = —A{fj). 

• Periodicity: A{fj ± T) = A{fj) . 

• Continuous intervals and singularities: A{fj — t- ±T/4) = ±oo . 

Most importantly, it is continuous only at a finite microscopically-small time interval 
T ~ 1/Aqcd and corresponds to spatially-homogeneous pulses of the gluon field potential 
with a constant energy-density. We notice also that once the compensation condition (13. 9p 
has been satisfied at a particular moment in time, it holds true for any later times, and the 
YM fields and large negative Ainst disappear from the resulting equations and do not partic- 
ipate in the Universe evolution any longer. We, therefore, arrive at quasistationary regime 
when the Universe evolution is completely determined by usual matter and uncompensated 
cosmological A-term only. 

What is the physical interpretation of the result (14.51) ? Obviously, such a solution with 
regular singularities does not have a quasiclassical interpretation (except for a vicinity of 
the midpoints of the continuous intervals where the gluon field potential is small and slowly 
changing). In practice, we deal with a sequence of quantum fiuctuations of the YM field in 
time, or, in fact, with the finite-time instantons. The creation and annihilation of such finite- 
time instantons should have essentially quantum nature like the QCD instantons. In order to 
regularize singularities in the quasiclassical solution (14. 5 p one needs to turn into a complete 
quantum theory taking into account quantum corrections due to e.g. fermion-antifermion 
pair creation and annihilation processes in the early Universe. 



V. CONCLUSION 



Based on the quasiclassical result (14. 5 P only, one can naively conjecture that the cosmolog- 
ical evolution of the YM field emerges a sequence of quantum tunneling transitions through 
the time harriers represented by the regular singularities in the solution (14. 5p . Thus, we have 
observed that the positive constant energy density of spatially-homogeneous finite-time in- 
stantons in the early Universe can be canceled with the negative constant QCD contribution 
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from spatially-inhomogeneous gluon field fluctuations induced by a similar quantum tunnel- 
ing of the gluon field, but through spatial (not time!) topological barriers between different 
classical vacua. This exhibits a remarkable similarity and interplay between instantons of 
different types in the early Universe evolution. 

From the point of view of quantum tunneling, the chain of quantum fluctuations in a 
certain approximation can be considered as non-linear oscillations. In order to regularize 
the infinite end-points of the continuous time intervals within the quasiclassical approach, 
one could therefore consider a continuous smearing of the resulting fluctuations by means 
of a non-linear continuous parameterization such as 

. , ^, , cos^fwr?) 
a sm[uiri) + b — — 

with adjustable parameters a, h and uj. This approximation has been qualitatively compared 
to the exact quasiclassical solution (14. 5 p and approaches it in the limit of small a — > 0, while 
the classical non-linear solution (12. 4 p is reached in the limit a — ?■ 6 (up to an arbitrary 
initial phase) with u being dependent on the initial amplitude. Thus, the parameterization 
(15. ip represents a simple continuous interpolation between the classical and quasiclassical 
solutions, and can be used in practical calculations in the quasiclassical limit of the theory. A 
detailed quantum analysis of the finite-time instantons is certainly a necessary step forward 
and is planned for a forthcoming study. 
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